Using Marie vacuum residue as raw material, the acetic acid⁃peroxyacetic acid oxidation system was used to pre⁃oxidize the vacuum residue, and N, N⁃dimethylformamide (DMF) was used to extract and remove sulfide in the residue. Under the conditions of m(residue)/m(solvent)=1 ∶ 1, Investigate the effects of oxidant dosage, reaction time and temperature on the removal of sulfur content in residue oil. X⁃ray photoelectron spectroscopy (XPS) and Fourier infrared spectrometer (FT⁃IR) The organic sulfur type distribution and functional groups of the residual oil before and after oxidation were detected. The results show that under the conditions of m(residue oil)/m(solvent)/m(acetic acid⁃peroxyacetic acid)=5 ∶ 5 ∶ 4, oxidation temperature 70 ℃, and oxidation time 50 minutes, the C-S stretching vibration absorption peak of the oxidized residue is obviously weakened. Thiophene type sulfur can be removed by extraction after being oxidized to sulfones. The sulfur content of the residue is reduced from 4.01% to 2.67%, the desulfurization rate reaches 33.37%, and the residue yield is 94.82%.

Through grey prediction GM (1, 1) model calculation and Matlab software prediction, the prediction of ground motion response of 400 m ³ storage tank is studied.First, the liquid height difference of the reservoir is set to 0.5 m, the total force, the base bending moment of the storage tank and the axial stress of tank wall under the horizontal earthquake are predicted under the different height of liquid storage.The prediction absolute error map is drawn, and the influence of earthquake intensity on the prediction error is obtained. Then, in order to further reduce the prediction error, the prediction effect is improved and the liquid height difference of the storage tank is adjusted to 1.0 m，reprediction the dynamic response of the storage tank. It is concluded that the interval of the high liquid difference data sequence of the storage tank increased from 0.5 m to 1.0 m. The relative error of the prediction results of the total horizontal seismic force and the base bending moment of the storage tank is reduced, after improvement, the prediction effect of bending moment and the axial stress of tank base are good.

The direct discontinuous Galerkin (DG) method was applied to calculate the soil instable temperature field of underground heated pipeline. This method applied partial integration from each computational cell for the heat transfer equation and obtained weak formulation.The direct DG was obtained by building appropriate numerical flux. The ordinary differential equations were gained by solving spatial discretization by means of the explicit Euler method. Through the calculation of the radial temperature field of heated pipeline ,its results coincide with other references'. Thus means that the direct discontinuous Galerkin method is an efficient and accurate numerical method for the calculation of instable temperature field.

The local discontinuous Galerkin method for the two-dimensional equation was introduced. By introducing the auxiliary variable, the second derivative heat conduction equation equations were rewritten into first-order partial differential equations, then the discontinuous Galerkin method was applied to the system. Finally the method yields ordinary differential equations systems. The explicit method was used for the time discretization. The numerical experiments were presented to verify the accuracy of the method. 